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Abstract 

Decoherence phenomena are pervasive in the arena of nanostructures but 
perhaps even more so in the study of the fundamentals of quantum mechan- 
ics and quantum computation. Since there has been little overlap between the 
studies in both arenas, this is an attempt to bridge the gap. Topics stressed 
include (a) wave packet spreading in a dissipative environment, a key element 
in all arenas, (b) the definition of a quantitative measure of decoherence, (c) 
the near zero and zero temperature limit, and (d) the key role played by initial 
conditions: system and environment entangled at all times so that one must 
use the density matrix (or Wigner distribution) for the complete system or 
initially decoupled system and environment so that use of a reduced density 
matrix or reduced Wigner distribution is feasible. Our approach utilizes gen- 
eralized quantum Langevin equations and Wigner distributions. 
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I. INTRODUCTION 



Quantum interference effects have become pervasise in two rather disparate areas of 
physics, (a) mesoscopic systems [1-9] and (b) fundamental quantum physics and quantum 
computation [10-24]. Whereas, of necessity, there is a certain commonality of approach, 
in general each area tends to emphasize the use of particular techniques. Since our own 
emphasis has been in the latter area, we feel that it might be worthwhile to present an 
overall account of our approach to an audience whose main interests lie in the former area. 

The generic problem of interest is the motion of a quantum particle in an environment. 
This problem has been addressed by a variety of methods such as the Feynman- Vernon 
functional integral approach, the use of master equations and various stochastic methods. 
We have found that the generalized quantum Langevin equation (GLE), supplemented by 
use of the Wigner distribution function, provides the basis of a powerful and physically 
transparent approach to all such problems. 

In Sec. II, we discuss fundamentals. Starting with the Hamiltonian for the whole system, 
we summarize how the GLE was derived and we discuss its properties. The solution for a 
non-localized quantum particle in a very general dissipative environment is then presented. 
Results are also presented for the correlation and commutation of both the coordinate and 
the fluctuation force operators appearing in the GLE. In Sec. Ill, we derive s(t), the mean- 
square displacement of the quantum particle in a dissipative environment [see (3.1) below]. 
The latter result, which is valid for arbitrary temperatures and arbitrary disspation, is a 
key ingredient of a very general expression for the spreading of the width of the particle's 
wave packet [see (3.11) below]. It also plays a vital role in quantum transport studies. In 
addition, it is relevant to a determination of the phase-breaking time r^, of interest for many 
mesoscopic systems. In Sec. IV, we apply our results to a calculation of the decoherence of 
interference terms in a quantum superposition. As a measure of decoherence, we define an 
attenuation coefficient a(t), the general result for which depends strongly on s(t) [as may be 
seen from (4.3) and (3.11) below]. Sec. V presents a generalization to the case of an external 
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field and, in Sec. VI, we present our conclusions. Fig. 1 presents a schematic summary of 
the progression from fundamental theory to applications. 



II. FUNDAMENTALS 

In recent years, there has been widespread interest in dissipative problems arising in 
a variety of areas in physics. As it turns out, solutions of many of these problems are 
encompassed by a generalization of Langevin's equation to encompass quantum, memory, 
and non-Markovian effects, as well as arbitrary temperature and the presence of an external 
potential V(x). As in [25], we refer to this as the generalized quantum Langevin equation 
(GLE) 

mx+ ! dt'n(t-t')x(t')+V'(x) = F(t) + f(t), (2.1) 

J — oo 

where V'(x) = dV(x)/dx is the negative of the time-independent external force and ji{t) is 
the so-called memory function. F(t) is the random (fluctuation or noise) operator force with 
mean {F(t)) = and f(t) is a c-number external force (due to an electric field, for instance). 
In addition, it should be strongly emphasized that "- the description is more general than 
the language -" [25] in that x(t) can be a generalized displacement operator (such as the 
phase difference of the superconducting wave function across a Josephson junction). 

A detailed discussion of (2.1) appears in [25]. In particular, it was pointed out that 
the GLE, while very general, can be realized with a simple and convenient model, viz., the 
independent-oscillator (10) model. The Hamiltonian of the 10 system is 

» = L + v{x) + £ (4 + W> ( « - *> 2 ) - < 2 - 2 > 
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Here m is the mass of the quantum particle while rrij and uij refer to the mass and frequency 
of heat-bath oscillator j. In addition, x and p are the coordinate and momentum operators 
for the quantum particle and qj and pj are the corresponding quantities for the heat-bath 
oscillators. The infinity of choices for the rrij and ojj give this model its great generality [25]. 
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Use of the Heisenberg equations of motion leads to the GLE (2.1) describing the time 
development of the particle motion, where 



ti(t) = ^2mjUjC08(ujt)6(t), (2.3) 

j 

is the memory function, with 6{t) the Heaviside step function. Also 

F(t) = J2m^(t), (2.4) 

j 

is a fluctuating operator force with mean (F(t)) = 0, where q h (t) denotes the general solution 
of the homogeneous equation for the heat-bath oscillators (corresponding to no interaction). 
These results were used to obtain the results for the (symmetric) autocorrelation and com- 
mutator of F(t), viz., 

l(F(t)F(t') + F(t')F(t)) 
1 f°° 

= -/ dcoRe[jl(uj + iO + )}hujcoth(hu/2kT)cos[u;(t-t')] (2.5) 
71" Jo 

\F(t), F(t')} = — / duRe{jl(u + iO + )}cu 
<~> Jo 

xsinu(t-t'). (2.6) 
Here fl(z) is the Fourier transform of the memory function: 

PCX) 

jl(z) = / dtfx(t)e izt . (2.7) 
Jo 

Equation (2.5) is often referred to as the second fluctuation-dissipation theorem and we 
note that it can be written down explicitly once the GLE is obtained. Also, its evaluation 
requires only knowledge of Rep,(u). On the other hand, the first fluctuation-dissipation 
theorem is an equation involving the autocorrelation of x(t) and its explicit evaluation 
requires a knowledge of the generalized susceptibility a(cu) (to be defined below) which is 
equivalent to knowing the solution to the GLE and also requires knowledge of both Re/2(u;) 
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and Im/t(u;). This solution is readily obtained when V(x) = 0, corresponding to the original 
Brownian motion problem [26]. As shown in [25], a solution is also possible in the case of 
an oscillator. Taking V(x) = \Kx 2 = \muj\x 1 , these authors obtained the solution of the 
Langevin equation (2.1) in the form 

x(t)= f dt'G(t-t'){F(t') + f(t')}, (2.8) 

J — oo 

where G(t), the Green function, is given by 

1 r°° 

G(t) = — dua(u + iO + )e- luJt , (2.9) 
with a(z) the familiar response function 



—mz z — izji{z) + K 
It is often convenient to write (2.8) in the form 



x(t) =x s (t)+x d , (2.11) 

where x d is the "driven" contribution [27] due to the external force f(t) and x s is the contri- 
bution due to the fluctuation force F(t). Here we have introduced a subscript s to emphasize 
that x s (t) is a stationary operator-process, in the sense that correlations, probability distri- 
butions, etc. for this dynamical variable are invariant under time-translation (t — > t + 1 ). 
In particular, the correlation, 

C (t - t') = ^(x s (t)x s (t') + x s (t')x s (t)) 

= - f dujlm{a(io + i0 + )}coth^- cos co(t-t'), (2.12) 
7r J 2kT 

is a function only of the time-difference t — t'. Furthermore, [27] 

C d (t,t') = ^{x(t)x(t')+x(t')x(t)) 

= C (t-t') + {x(t)){x(t')), (2.13) 



where (x(t)) is the steady mean of the driven motion. 
Also, taking the Fourier transform of (2.8), we obtain 



a(uj){F(u) + f(u)}, 



(2.14) 



where the superposed tilde is used to denote the Fourier transform. Thus, x{uS) is the Fourier 
transform of the operator x(t): 



We have now all the tools we need to calculate observable quantities. It is also useful to 
note that the commutator, which is temperature independent, is given by the formula [28] 



III. MEAN SQUARE DISPLACEMENT AND WAVE PACKET SPREADING 

The mean square displacement of the quantum particle in a dissipative environment, 
s(t) say, plays a key role in all of our subsequent discussions. In particular, it determines 
the diffusion time through a mesoscopic system of interest. In order to concentrate on the 
fundamental features, we take V = and f(t) = but we will return to f(t) ^ in Section 
5. Then using (2.12), we obtain 




(2.15) 




(2.16) 



S (t) EE Qx B (t) - X S (0)] 2 > 



2{C (0)-C (t)} 




(3.1) 



where now 
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(3.2) 



mz 2 — izp,(z) 
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It is informative to consider some limiting cases. In particular, in the quantum 
physics/quantum optics literature the Ohmic model (no memory in the GLE), for which 



jl(ui) = rwy = constant (Ohmic), (3.3) 

is often considered. Then, in the high temperature limit one obtains [29], with t > 0, 



2kT 

s(t) = {e-T* - 1 + t*} , kT » fry. (3.4) 



For long times 



2kT 

s(t) -> 1, 7 t » 1 (3.5) 

which is the familiar Einstein relation from Brownian motion theory. However, for short 
times (which are more characteristic of decoherence decay times [16]) 



kT 

s(t) -> — t 2 , -ft « 1, (3.6) 
m 

independent of 7. Thus, in this simple case, the familiar diffusion coefficient 



1 kT 

D=- 2 m = -, (37) 

so that the phase-breaking time = 7 _1 is given by 

r* = (3.8) 

More generally, it is clear from the above that is a complicated function of s(t), depending 
on the temperature, on the nature of the environment and on the time scale of interest. 

In particular, the situation is very different for low temperatures (kT << fry), in which 
case the main contribution is from the zero-point (T = 0) oscillations of the electromagnetic 
field. This calculation has recently been carried out [24] and it was found necessary to use a 
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model incorporating memory since the Ohmic model lead to singular and unphysical results 
for very short times. In particular, the single relaxation time (r) model 



1 — zt j 

leads to the result [24] 



m = (3.9) 



s(t) = -^-t 2 (log^ + lE - \ ) , T«t« 7"\ 
Txm i m 2 ) 



(3.10) 



where 7# is Euler's constant. It would be of interest to obtain corresponding results for the 
actual dissipative models which describe various mesoscopic systems. 

Whereas s(t) is of interest per se, it also plays a key role in the determination of wave 
packet spread and attenuation of coherence. Recently [16], we obtained a general formula 
for the spreading of a free particle Gaussian wave packet: if initially the complete system is 
in equilibrium at temperature T and the wave packet describing the quantum particle has 
a width a, then after a time t the width is given by 

2 _ [x(0),x (f ff + 

where [x(0), x(t)] is the commutator. For example, for a free particle in a high temperature 
environment 

^ 2 W^ 2 + 7^ + — t 2 . (3.12) 

For T = 0, (3.12) reduces to the familiar formula of quantum mechanics [30,31]. In general, 
the width w(t) is a vital ingredient in the determination of decoherence decay times, as we 
shall see in the following section. 
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IV. DECOHERENCE 



Decoherence refers to the destruction of a quantum interference pattern and is relevant 
to the many experiments that depend on achieving and maintaining entangled states. Ex- 
amples of such efforts are in the areas of quantum teleportation, [13] quantum information 
and computation, [10,11] entangled states, [12] Schrodinger cats, and the quantum-classical 
interface. For an overview of many of the interesting experiments involving decoherence, we 
refer to [12-15]. 

Much of the discussion of decoherence [14,16-24] has been in terms of a particle moving 
in one dimension that is placed in an initial superposition state (Schrodinger "cat" state) 
corresponding to two widely separated Gaussian wave packets, each with variance a 2 and 
separated by a distance d. For such a state the probability distribution at time t can be 
shown to be of the form: 

P(X ' t] = 2(1 + e-*W) { P °0 " 5'*) + P ° (X + P> 

+2 e-/-(.» oW P o( ,, ( )oosM^}, (4.1) 

where Po is the probability distribution for a single wave packet, given by 

p ^^7^m eM ^ 1 (4 ' 2) 

Here and in (4.1) w 2 (t) is the variance of a single wave packet, which in general is given 
by (3.11). Also a 2 is the initial variance and s(t) is given by (3.1). In (4.1) the first two 
terms within the braces correspond to the two wave packets, centered at ±d/2, expanding 
independently, while the third term is the interference term. Decoherence refers to the 
destruction of interference, a measure of which is given by the attenuation coefficient a(t) 
which can be defined as the ratio of the factor multiplying the cosine in the interference 
term to twice the geometric mean of the first two terms [16,17] 

In general, the interference term disappears in time due to either temperature or dissi- 
pative effects. There are various measures of decoherence, based on decay of diagonal and 
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off-diagonal density matrix elements or probability distributions in phase space, momentum 
space or coordinate space [20] but we consider the latter to be the most desirable because 
it is closest to experiment. Thus, we measure the disappearance of the interference term, 
that is, the loss of coherence (decoherence), by calculating a(t), leading to the very general 
result [16]. 



a(t) = exp 



s(t) 



(4.3) 



8a 2 w 2 (t) 

As a concrete example, if we consider a low temperature environment then, using (3.10), 
(3.11) and (4.3), and we obtain [24] 



a(t) = exp 



log — + 7£ - - 

m 2 



with 



r « t « (C/m) 1 



(4.4) 



ma 



18^ 



TO = 



(4.5) 



d V K 

which leads to very short decoherence times in the case of the macroscopic separations d 
generally considered in the quantum physics/quantum optics arena. However, for nanostruc- 
tures, the latter assumption is no longer valid (simply because values of d are microscopic) 
and, concomitantly, the short-time approximation ('-ft << 1) will no longer be valid in many 
cases simply because values may be comparable to r^. 



V. PRESENCE OF AN EXTERNAL FIELD (f(t) + 0) 

As we saw in section 2, in particular (2.11) and (2.13), the effect of an external classical 
force f(t) is to replace x s — > x s + x d and C — > Cd- This leads, with the help of (3.1) to the 
conclusion that, in the expression for a(t) given by (4.4), we simply replace s(t) — > s(t) + s d , 
where 
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Sd(t) = (x 2 d (t)) 

t ft 

dt' dt"G(t-t')G(t-t")g(t' -t"), (5.1) 
Jo 

with 

g(t>-t") = (f(t>)f(t")). (5.2) 

In order to proceed further, one must, of course, specify the time dependence of the external 
force. For details, we refer to [23]. 

VI. CONCLUSION 

The material present above is an overview of our work in the arena of quantum 
physics/quantum optics, where the models we have discussed display great generality: ar- 
bitrary temperature and arbitrary heat bath allowing for non-Markovian (memory) effects 
and with no restriction to initial decoupling of the bath and system. They are relevant not 
only to fundamental questions but also have implications for quantum computation studies 
(since decoherence is always a stumbling block to the maintenance of superposition states). 
Our hope is that they will stimulate new ideas in the study of nanostructures and mesoscopic 
systems generally since analytic results have been obtained in most cases which facilitate 
attainment of physical insight into the important facets of the various problems. The next 
steps should include, in particular, an analysis of more complicated potentials and also heat 
baths which display memory effects obtained from physical models. The latter, in turn, 
will have to reflect, in many cases, the actual systems of interest since, for example, the 
inelastic phase coherence time has very different forms depending on the parameters of the 
system [8,9]. Furthermore as mentioned previously, consideration of wave packets separated 
by nano distances will necessitate the consideration of a greater range of decoherence decay 
times with the consequence that many assumptions inherent in existing calculations will no 
longer be valid. 
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